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ABSTRACT 

We present 2 cm and 3.6 cm wavelength very long baseline interferometry images of the compact 
radio continuum sources in the nearby ultra-luminous infrared galaxy Arp220. Based on their radio 
spectra and variability properties, we confirm these sources to be a mixture of supernovae (SNe) and 
supernova remnants (SNRs). Of the 17 detected sources we resolve 7 at both wavelengths. The SNe 
generally only have upper size limits. In contrast all the SNRs are resolved with diameters > 0.27 pc. 
This size limit is consistent with them having just entered their Sedov phase while embedded in 
an interstellar medium (ISM) of density 10 4 cm -3 . These objects lie on the diameter-luminosity 
correlation for SNRs (and so also on the diameter-surface brightness relation) and extend these 
correlations to very small sources. The data are consistent with the relation L oc D~ 9 / 4 . Revised 
cquipartition arguments adjusted to a magnetic field to relativistic particle energy density ratio of 1% 
combined with a reasonable synchrotron-emitting volume filling factor of 10% give estimated magnetic 
field strengths in the SNR shells of ~ 15-50 mG. The SNR shell magnetic fields are unlikely to come 
from compression of ambient ISM fields and must instead be internally generated. We set an upper 
limit of 7 mG for the ISM magnetic field. The estimated energy in relativistic particles, 2%-20% of 
the explosion kinetic energy, is consistent with estimates from models that fit the IR-radio correlation 
in compact starburst galaxies. 

Subject headings: galaxies: individual (Arp220) — galaxies: starburst — ISM: supernova remnants 



1. INTRODUCTION 

Observations of radio supernovae (SNe) and super- 
nova remnants (SNRs) provide an important means to 
study astrophysical processes occurring in dense nuclear 
starbursts. Radio observations of these objects are free 
from the dust obscuration that hamper observations at 
shorter wavelengths while the high angular resolution af- 
forded by very long baseline interferometry (VLBI) ob- 
servations cannot be matched by any other technique. 
Monitoring the rate of appearance of radio SNe can po- 
tentially constrain the stellar initial mass function (IMF) 
and check whether it is modified in e xtremely dense en - 
vironments such as found in Arp220 (|Parra et alj [2007). 
The SNRs are the acceleration sites of the relativis- 
tic particles that give rise to radio emission in star- 
forming regions and hence their properties are central 
to un derstanding the FIR-radio correlation ([Lacki et al.l 
2010). Finally, radio SNe and SNRs can be used as in 
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situ probes to constrain the interstellar medium (ISM) 
properties such as density, pressure, and magnetic field 
strength. Recent papers dealing with VLB! observations 
of SN e/ SNRs in starbur s t gala xi es include iFenech et al.1 
&il(h. ILenc fc Tingavl (l2009h. lUlvestadl (f2009fl . and 
iPerez-Torres et al.l (|200a ). 

Arp220 is the nearest and best studied ultra-luminous 
infrared galaxy (ULIRG). For over a decade it has been 
the subj ect of a glo b al VL BI campaign at cm wave- 
lengths. Smit h et al.l (|1998f ) made the first detection of 
a number of compact sub-parsec sized sources at 18 cm. 
These sources have been monitored and new objects dis- 



Rovilos et afl 12005; 



2007) reported the 



covered in successive 1 8 cm epochs 
lLonsdale et "all 120061 ). iParra et al.l 
first detections of these objects at the shorter wave- 
lengths of 13 cm, 6 cm, and 3. 6 cm. Based on th e source 
radio light-curves and spectra, IParra et"aTI (|2007t ) argued 
that the compact radio sources in Arp220 comprise a 
mixed population of radio SNe and SNRs. 
This paper presents the results and analysis of new 
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TABLE 1 

High Frequency VLBI Observations of Arp220 



Epoch 


Code 


A 


Array 


ff rms 


Beam Size 






(cm) 




(/ijy beam -1 ) 


(mas 2 ) 


2006.02 


BP129 


13.26 


VLBA 


129.54 


3.6x6.6 






6.02 a 




41.43 


1.8x3.5 






3.56 




86.73 


1.7x3.1 


2006.91 


GC028A 


3.56 


Global VLBI 


34.7 


0.56x1.34 


2006.99 


GC028B 


2 


HSA 


28.12 


0.43x1.26 


2008.44 


GC031A 


6.02 


Global VLBI 


12.58 


0.72x2.1 



Note. — a Rc-rcduccd 6 cm observations from experiment BP129. 



high sensitivity and resolution VLBI observations at 
3.6 cm and 2 cm wavelengths. The main objectives of 
these new observations were to extend source spectra to 
higher frequency, look for high-frequency variability, and 
most importantly attempt to spatially resolve sources or 
set limits on source sizes. In Section [2] we present details 
of the high-frequency observations and their reduction; 
we also describe other ancillary VLBI data used in our 
analysis. In Section [3] we describe our results including 
the source detection and resolution criteria, the estima- 
tion of source sizes and of source spectra. In Section|4]we 
discuss our results and their implications for the physics 
of the sources and the ambient ISM. Finally in Section [5] 
we give our conclusions. In this paper we assume a dis- 
tance of 77 Mpc to Arp220 at which an angular size of 
1 mas corresponds to 0.37 pc. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. 3.6 cm Global VLBI 

Arp220 was observed at 3.6 cm (8.4 GHz) on 2006 
November 28 as part of European VLBI Network (EVN) 
+ Very Long Baseline Array (VLBA) global VLBI exper- 
iment GC028A. The experiment used the 10 stations of 
the VLBA, the Green Bank Telescope (Gb), the phased 
Very Large Array (VLA; Y27), Arecibo (Ar), and five 
EVN stations (Ef, On, Mc, Nt, Wb). The observations 
were performed using a dual polarization 256 Mbit s _1 
recording mode. The data were correlated at JIVE in the 
Netherlands in a single pass using a phase center between 
the two nuclei and high spectral and time resolution to 
allow wide-field imaging. Unfortunately, no fringes were 
detected at the Green Bank Telescope which significantly 
reduced the sensitivity of the observations. Also due to 
a scheduling error at the station Ar was only available 
for a very small fraction of the time. 

The data were reduced using AIPS in a standard man- 
ner. Initial amplitude calibration was carried out using 
monitored system temperatures and station gains. We 
estimate that absolute amplitude calibration should be 
accurate to ^5%. The main phase calibrator was the 
nearby (0?5 separation) source J1532+2344 which was 
observed in a rapid switching cycle with Arp220. Simul- 
taneous CLEAN images of size 4096 pixels x 2048 pixels 
(pixel spacing 0.1 mas) were then made of the eastern 
and western nuclei using robust data weighting to give a 
good compromise between sensitivity and resolution. Us- 
ing the initial phase-referencing phase calibration it was 
found impossible to make thermal noise limited images. 
To improve the dynamic range we performed phase self- 
calibration on the Arp220 data itself using long solution 



times (« 30 minutes) to obtain sufficient signal to noise. 
The resulting slowly varying phase solutions were con- 
sistent with unmodeled atmosphere delays differentially 
affecting target and calibrator. Finally CLEAN was ap- 
plied again and noise limited images produced (see Ta- 
ble [TJ . The final images are shown in the central panels 
of Figures [T] and [2j with inset panels showing detailed 
images of individual sources. In the rest of this paper we 
refer to experiment GC028A when mentioning the 3.6 cm 
data or image if not otherwise specified. 

2.2. 2 cm HSA 

Arp220 was observed at 2 cm (15 GHz) in a full track 
observation on 2006 December 28 using the High Sen- 
sitivity Array (HSA; expt GC028B). The experiment 
used the 10 stations of the VLBA, the Green Bank Tele- 
scope (Gb), the phased VLA (Y27), and the Effelsberg 
(Eb) antenna. The observations were performed in a 
256 Mbit s _1 dual polarization mode. The data were 
correlated in Socorro, New Mexico in one pass with a sin- 
gle phase center located midway between the two nuclei 
and high spectral and time resolution to allow wide-field 
imaging. The data were reduced using AIPS in a stan- 
dard way similar to that described in Section 12.11 After 
experimentation with many uv data weighting schemes 
pure natural weighting of the visibility data was cho- 
sen. This choice gave the highest possible sensitivity at 
the expense of resolution; however for other choices most 
sources were not detectable. As for our 3.6 cm wave- 
length data, phase self-calibration of the Arp220 data 
was required to achieve a thermal noise limited image 
(see Table [l}. Our resulting 2 cm wavelength images for 
individual sources in each nucleus are shown in insets in 
Figures Q] and [2j respectively. 

2.3. Other VLBI Data 

A number of ancillary VLBI data sets, providing mul- 
tifrequency monitoring information, were also used in 
the analysis presented in this paper. These observations 
were used to primarily classify sources as SNe or SNRs 
(see Section l4"7Tj) . A full presentation of these monitoring 
data, including data on many sources only detected at 
lower frequency, will be given in a future paper (Batcjat 
et al. 2012, in preparation). 

The most extensive monitoring data exist at 18 cm 
wavelength where a total of nine epochs spanning from 
1994.87 to 2006.43 have been observed and reduced. 
These data provide information on source long wave- 
length variability on decade timescales. We have re- 
analyzed the complete 18 cm data set in a consistent 
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fashion carefully correcting for epoch dependent varia- 
tions in flux scale. Additionally to this, shorter wave- 
length variability can be searched for by comparing the 
res ults at 6 cm from epoch 2006.02 (experiment BP129, 
see lParra et al.l2007t ) re-reduced by us in order to achieve 
a lower er rms (see Table [2] for the new flux densities) with 
data 2.42 years later at epoch 2008.44 from a recently 
reduced observation (experiment GC031A, see Table [J). 
This latter epoch is one of a series of 6 cm monitoring 
experiments which we are presently reducing. Finally we 
compare ou r new GC028A 3.6 cm data with data from 
IParra et al.l (|2007T) taken 0.9 years earlier. Because of the 
low SNR of this old 3.6 cm wavelength data it has not 
been possible to re-reduce it using self-calibration tech- 
niques (unlike 6 cm data from the same experiment). 

3. RESULTS AND ANALYSIS 

3.1. Source Detection 

Visual inspection of the final 2 cm and 3.6 cm images 
(Figures Q] and [2]) shows many sources clearly visible at 
both wavelengths. To more rigorously define a list of de- 
tected s ources we applied a two-pass search method as 
used by IParra et"aTI J2007). First the whole east and 
west nucleus images were searched for sources signifi- 
cantly above the noise, using a detection limit which min- 
imized the chance of false detections. Then in a second 
pass small regions were searched around the positions of 
known sources (previously detected at 18 cm, 6 cm and 
3.6 cm since the f irst detection of compact sources by 
ISmith et all lfl998|n using a lower detection threshold. 
For a detection limit set to rja where a is the noise rms, 
the probability of one or more false detections in an im- 
age is F = 1 — P(I < rja) Ns , where P(I < rja) is the 
cumulative probability of a noise spike at a given beam 
area being less than r/a and N s is the number of beam 
areas searched. Histograms of the pixel values in source 
free areas at both wavelengths were close to Gaussian 
out to the most extreme pixel values, hence P could be 
calculated assuming Gaussian statistics. In both passes 
we chose rj such that F was 0.2%. Given that we have 
two passes and two wavelengths, this choice gave a final 
probability of < 1% of getting one or more false detec- 
tions. 

For the first pass, given the number of beam areas 
at 2 cm, the above chosen value of F corresponded to 
?7 = 5.6. Although the number of searched beams was 
slightly less at 3.6 cm wavelength we conservatively used 
the same detection criteria. This first pass applied sep- 
arately to both wavelengths resulted in the detection of 
14 sources in the western nucleus (two of which, namely 
W58 and W60, were new detections) and 3 in the east- 
ern nucleus. All of these sources were detected at both 
wavelengths. In the second pass the search regions were 
limited to boxes of 4 mas 2 centered on the positions of 
known sources. Given the smaller number of beam areas 
searched in this pass the critical detection threshold was 
smaller (r) = 3.8). Despite this lower limit no additional 
detections were made. 

Table [2] gives the absolute positions of the 17 detected 
sources, estimated by fitting a Gaussian to the 2 cm im- 
age of each source. Given the relatively close-by cali- 
brator source used (0?5 dista nt) the absolute as tromet- 
ric accuracy is likely limited (jPradel et al.l 120061) by the 



accuracy of the phase cal ibrator position (e stimated for 
J1532+234 to be 0.5 mas; lPetrov et aH l2005l). In Tabled 
we therefore give positions rounded to the nearest mil- 
liarcsecond in each coordinate. Columns 9-12 of Table[2] 
give source peak and total flux densities at 3.6 cm and 
2 cm from experiments GC028A and B. The total flux 
densities were measured within a tight box surrounding 
each source. For other epochs and wavelengths (Columns 
5-8) only peak flux densities are given; since however 
sources are unresolved in these data, peak brightness 
measurements should also be good estimates of source to- 
tal flux densities. Note that the 6 cm flux density values 
given for BP 129 in Tableware based on a re-reduction of 
the data using self-calibrati on while the value s at other 
wavelengths are taken from IParra et al.l (|2007l) . For the 
re-reduced 6 cm data source flux den sities are on averag e 
a factor of 1.36 larger than found by IParra et al.1 ((2007); 
this increase is likely due to an increase in phase coher- 
ence of that data. A possible slight overestimate in flux 
densities due to noise biasing of low SNR self-calibrated 
data potentially could also be present but this potential 
effect is hard to quantify without detailed simulations. 

3.2. Size Estimation 

A primary goal of our new high-frequency observations 
was to measure or set limits on the source sizes. From 
visual inspection of the inset images in Figures [T] and [2j 
there are several candidate resolved sources (i.e., W18, 
W33, W42, and E14) where the 50% of peak contour, 
at least at one wavelength, encloses a significantly larger 
area than the beam and so appear resolved. Quantita- 
tive tests are however required to confirm or reject these 
visual impressions and to give size estimates with confi- 
dence limits. Two methods were developed to do this. In 
the method described in Section 13.2.1 1 we tested whether 
the null hypothesis that a source was unresolved could be 
rejected while in Section 13.2.21 we describe a procedure 
applied to give best estimates (with error bars) of each 
source's size. In Section 13.2.31 we summarize the results 
obtained after applying these methods to our data. 

3.2.1. Test of Resolution 

We tested whether or not each source was consistent 
with the null hypothesis of a point source convolved with 
the CLEAN Gaussian restoring beam (having minor and 
major axis dispersions of <r x and a y ). For each source we 
calculated an observed minor axis dispersion (T XobB from 

2 J J I(x,y)Ax 2 dxdy 
CTx ° bs J J I(x,y)dxdy ' 1 ' 

where Aa; is measured along the minor axis through the 
source center. If, given the expected noise, a XobB was 
significantly larger than er x the point-source hypothesis 
was rejected. Analytic calculation of how much the for- 
mer quantity must exceed the latter to be confident of 
resolution is difficult, first because of pixel-to-pixel corre- 
lations in noise, and second because of blanking applied 
to the data. Any practical method of estimating <r Xobs 
must first blank the image I(x, y) below say, 3tJ rms or 
else the statistic is dominated by distant noise peaks. 
This blanking however makes the problem non-linear. 

To circumvent the above problems we applied a statis- 
tical bootstrapping technique. For each detected source a 
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Fig. 1. — Central panel shows the 3.6 cm image of the wes tern nucleus of Arp220 spanning a region 410 mas by 205 mas. Fourteen sources 
were detected above a 5.6cr detection threshold (see Section |3. Il l at both 2 cm and 3.6 cm. The surrounding panels show zoomed images of 
these sources at 2 cm (left panels) and at 3.6 cm (right panels). The same color bar is used for the central panel and all zoomed images at 
both 2 cm and 3.6 cm. Each zoomed image covers a region of 2.1 mas by 4.1 mas. The source number is displayed in the top left corner of 
the 2 cm zoomed image. Where the source number is displayed in red this means that the source is resolved at both 2 cm and 3.6 cm. A 
green number means resolved at 2 cm only, a cyan number means resolved at 3.6 cm only while a white number means unresolved at both 
2 cm and 3.6 cm. For more information about how resolution was determined see Section l3.2.1l The zoomed images are blanked at the 3cr 
level and contour levels are shown at 50%, 75%, and 95% of the peak intensity of each source (given in Table[2]l. The 50% contour of the 
CLEAN restoring beam (see Table [TJ at each wavelength is shown in the bottom left zoomed image. 



point of the same flux density convolved with the restor- 
ing beam was added at multiple (source- free) positions on 
the final image to test the effect of different realizations of 
the noise. Each realization was blanked below 3cr rms and 
the minor axis spatial dispersion <7 Xobs measured. The 
resulting histogram of measured dispersion values was 
then compared to that measured for the source using the 
same blanking criteria. If this latter quantity was greater 
than 90% of the histogram values the null hypothesis of 
a point source was provisionally rejected. 

Table [3] lists for each source and wavelength the per- 
centage probability of the point source hypothesis being 
rejected. According to our adopted criteria of 90% con- 
fidence 13 sources are provisionally resolved at 2 cm and 
8 at 3.6 cm. In total seven sources fulfill our resolution 
criterion at both bands. Since all these have < 1% cumu- 
lative probability of achieving their measured dispersion 
by chance if a point source we argue that for these sources 
the point source hypothesis can be rejected. 

3.2.2. Estimating Source Sizes 

Best estimates of source sizes with error bars were 
made using a variation on the method presented in Sec- 
tion 13.2.11 Again we added a model source to the final 
image at multiple positions to take into account different 



realizations of the noise. In this case however the model 
added was a spherical shell convolved with the CLEAN 
restoring beam, whose outer diameter we then adjusted 
at each position till the measured dispersion after blank- 
ing equaled the observed source dispersion after blank- 
ing. At positions where even a source of zero size was un- 
able to produce a minor axis dispersion small enough to 
equal that measured from the data we set the estimated 
source size to zero. The SNR shell model used had a ratio 
of inner to out er diameter 0.8 similar t o that measured 
for SN1993J (iMarti'-Vidal et al.l l201lh and SN2008iz 
( Brunthal er et al.l I2010T ) and cons istent with the obser- 
vation of evolved SNRs in M82 (TKronberg et all 119851: 



-onberg 
20061 : II 



Fenech et al.1 



Muxlow et al.1 119941 : iBeswick et all ._ 

20qH2010). 

From the resulting histogram of outer shell 

diameters we then obtained a best estimate of the source 
size and estimated error bars. In some cases we obtained 
good fits for the outer diameter (> 0) in over 90% of 
realizations. In these cases we took the best estimate of 
source size as the mean of the histogram over all real- 
izations with failed solutions set to zero size and defined 
error bars from its 10th and 90th percentiles. In other 
cases with < 90% of realizations giving good solutions 
we still set a size upper limit at the 90th percentile of 
fitted sizes but we set the lower limit to zero. 
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Fig. 2. — Central panel shows the 3.6 cm image of the ea ster n nucleus of Arp220 spanning a region 410 mas by 205 mas. Three sources 
were detected above a 5.6cr detection threshold (see Section |3.1|I at both 2 cm and 3.6 cm. For information on labeling, contours, and other 
image properties see Figure [T] 



TABLE 2 

Position of Detected Radio Sources and Flux Densities 







O2000 


<52000 


BP129 


BP129 


GC031A 


BP129 


GC028A 


GC028A 


GC028B 


GC028B 


Name 


SN 


15 h 34 m ... 


23°30'... 


13 cm 


6 cm 


6 cm 


3.6 cm 


3.6 cm 


3.6 cm 


2 cm 


2 cm 










Peak 


Peak 


Peak 


Peak 


Peak 


Integrated 


Peak 


Integrated 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Error (it) 








±130 


±41 


±13 


±87 


±35 


±35 


±28 


±28 


Wll 




57! 2299 


1T'502 


573 


581 


818 


357 


289 


279 


190 


163 


W12 




57? 2295 


11^524 


1011 


711 


953 


311 


389 


416 


258 


253 


W15 




57? 2253 


ll!'483 


579 


959 


1115 


707 


546 


837 


314 


500 


W17 


6 


57? 2241 


ll!'520 


477 


524 


454 


484 


218 


385 


150 


131 


W18 


7 


57? 2240 


H^'547 


559 


834 


738 


151 


272 


560 


156 


288 


W25 




57? 2222 


lK'501 


1069 


1141 


1479 


648 


752 


965 


390 


324 


W33 


11 


57? 2200 


11^491 


582 


417 


321 


397 


310 


409 


162 


182 


W34 




57?2195 


ll!'492 


699 


1066 


1458 


743 


650 


758 


263 


358 


W39 


12 


57?2171 


lK'485 


749 


428 


466 


120 


190 


247 


129 


91 


W42 


13 


57?2122 


ll!'482 


743 


693 


706 


574* 


273 


492 


240 


251 


W55 




57? 2227 


ll!'482 


121 


1000 


569 


1147 


822 


1001 


595 


693 


W56 




57? 2205 


11^491 




905 


1434 


750 


657 


815 


426 


446 


W58 




57? 2194 


lK'508 




408 


723 




322 


296 


226 


250 


W60 




57?2276 


ll!'546 




285 


1016 




310 


268 


200 


128 


E10 




57? 2915 


ll!'335 


227 


1034 


1631 


988 


1109 


1394 


577 


596 


Ell 




57? 2913 


ll!'333 




429 


310 




201 


171 


165 


195 


E14 




57 s 2868 


IE" 298 




698 


608 


549 


343 


509 


264 


400 


Note. — 


Col. (1) 


source name 


from [Lonsdale et al. 


(2006) for all sources 


except W55 and W56, which arc 


from IParra et aTJ (1200711 . 



Sources W58 and W60 are newly d etected so urces. W indicates t he source is located in the western nucleus while E stands for east. 
Col. (2): Name used in lSmith et all i 1998ft and lRovilos et alj d2005ft . Cols. (3) and (4): J2000 Right Ascencion and Declination obtained 
by fitting a gaussian to the sources in the highest frequency map (2 cm map). Col. (5) : 13 cm peak flux density (in fijy beam -1 ) 
from the observations performed in experiment BP129 and presented in Parra et al. (2007). Col. (6): re-reduced 6 cm peak flux density 
(in fj,Jy beam -1 ) from experiment BP129. Re-working of these data allowed us to produce a deeper map of both nuclei with cr rms — 
41.43 juJy beam -1 and to detect two more sources in the western nucleus (W58 and W60) and one more source in the eastern nucleus 
(Ell). Col. (7): Newly acquired (epoch 2008.44) 6 cm peak flux density (in fijy beam -1 ) which will be presented and discussed in 
detail in Batejat et al. (2012) (in pre paration). Col. (8) : 3.6 cm peak flux density (in fijy beam -1 ) from the observations performed 
in experiment BP129 and presented inParra et al. (2007). For source W42 (marked by an asterisk) we give the integrated flux density 
because this source appeared resolved. Cols. (9) and (10): 3.6 cm peak and integrated flux densities (respectively in fijy beam -1 and 
in /^Jy) from experiment GC028A discussed in this paper. Cols. (11) and (12): 2 cm peak and integrated flux densities (respectively in 
[iJy beam -1 and in //Jy) from experiment GC028B discussed in this paper. 
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3.2.3. Summary of Results for Source Resolution/Size 

The results of the statistical tests described in Sec- 
tions |32j] an d GL22] are summarized in Figure Three 
sources (W12, W25 and W60) are very compact and un- 
resolved at both wavelengths, using both resolution/size 
estimation methods. This demonstrates that residual at- 
mospheric phase errors do not have any appreciable effect 
in broadening our sources, since such errors would affect 
all sources. Seven sources (E10, Wlf , W34, W39, W55, 
W56, and W58) are provisionally resolved at one wave- 
length (usually 2 cm) and a further seven sources (Ell, 
E14, W15, W17, W18, W33, and W42) are resolved at 
both wavelengths. All of the latter group have diameters 
greater than 0.7 mas. 

Figure [3] shows in general a strong correlation between 
the sizes/limits measured at the two bands. Two clear 
exceptions however are W34 and W58, both of which are 
well resolved at 2 cm but are unresolved at 3.6 cm. In- 
spection of the detailed image for W34 at 2 cm shows 
that an isolated probable noise feature to the north-west 
gets through the initial 3cr blanking and contributes to 
a large minor axis spatial dispersion. The case of W58 
is less clear since there is no distinct peak off the main 
source just a broadened source with position angle differ- 
ent to the CLEAN beam. This could however be caused 
by a large noise peak that lies very close or on top of 
the source. The presence of one or two such anomalous 
sources is not unexpected by chance given the blanking 
method used. There are 17 sources each at two frequen- 
cies with approximately 10 independent beam areas per 
source box in Figures [1] and Given this total area 
there is a probability of 37% of detecting one or more 
such > 3cr noise peaks and a 8% chance of detecting two 
or more such peaks. 

It should be noted that although well resolved at both 
bands both W42 and W33 seem to be significantly larger 
at 3.6 cm than 2 cm; possible explanations for this are 
discussed in Section l4~2l In addition the estimated size of 
W42 at 3.6 cm in GC028A (ring outer diameter 0.82 mas) 
is much less than t hat claimed at the same frequency by 
IParra et all (|2007f ) from earlier BP 129 VLB A observa- 
tions (equivalent to an outer shell diameter of 3.6 mas 
after converting from the fitted Gaussian FWHM). It 
should be noted that the beam area of these earlier 
VLBA observations was much larger, by a factor of six, 
compared to our new global 3.6 cm observations. A con- 
sequence is that a shell of the size and brightness claimed 
from the BP 129 observations would be below the ther- 
mal noise per beam of full resolution GC028A data and 
so would be undetectable. To check this possibility lower 
resolution versions of the GC028A images and recent 
very sensitive global 6 cm observations (GC031A) were 
inspected but neither were consistent with the large shell 
structure claimed from the BP129 3.6 cm data. 

The ini tial argument for the resolution of W42 at 
3.6 cm in IParra et al.l (|2007l) was based on having two 
large VLBA beams areas above 50% of the peak source 
brightness. Significant residual atmospheric phase er- 
rors in the data (expected to be larger than those known 
to be present at 6 cm) may have contributed to an ap- 
parent source broadening, unfortunately the SNR is too 
low at 3.6 cm to allow self-calibration reprocessing of 
the BP 129 to quantify atmospheric broadening at this 



wavelength. Residual phase errors at 3.6 cm would have 
broadened all sources equally but when combined with 
the low surface brightness of W42 this might explain the 
BP 129 result. From these data the 50% of peak contour 
level is only thr ee times the rms of the thermal noise 
(see Figure 4 of IParra et al.1 (120071) ). Given t here were 
a total of 18 sources detected by IParra et al.1 (|2007l ) the 
probability of at least one of them having a 3a noise 
peak adjacent to the source, so giving the impression 
of a highly resolved source, is not negligible (2.5%) and 
furthermore this probability of false resolution would in- 
crease rapidly in the presence of even a moderate amount 
of atmospheric source broadening. 

3.3. Source Spectra 

The short wavelength observations presented in this 
paper allow us to extend the s ource radio spectra 
first studied by IParra et"al~l (|2007h to higher frequency 
(15 GHz, 2 cm). In addition the 3.6 cm observations and 
new 6 cm observations (see Table [lj can be compared 
to those in IParra et al.l (|2007|) to look for high frequency 
variability. Both spectral shapes and variability are use- 
ful diagnostics when attempting to classify sources as 
SNRs or SNe (see the Appendix) . Figure [4] shows the 
source spectra with the points measured in the period 
2006.02-2006.43 plotted in blue and those observed in 
the period 2006.91-2008.44 plotted in red. The measured 
flux densities at each epoch are taken from Tabled 

4. DISCUSSION 

4.1. Source Properties and Identification 

IParra et al.l ()2007f ) argued that the compact radio 
sources in Arp220 comprise a mixed population of SNe 
and SNRs; the former embedded in an ionized circum- 
stellar bubble and the latter strongly interacting with the 
surrounding ISM. In addition to the above two classes 
SN/SNR "transition" objects may also exist in which 
ISM interaction has begun but the swept up mass is less 
than the ejecta mass. Other compact sources might be 
associated with active galactic nuclei (AGNs). Below 
we discuss the expected properties and likely members 
of each class. The Appendix gives a description of the 
variability and spectral properties for each source and its 
classification. 

Supernovae. In these objects the synchrotron-emitting 
blast wave transits the dense circumstellar medium of 
the progenitor star, originating from a pre-explosion stel- 
lar wind with a r~ 2 density profile, which is ionized by 
the SN explosion. As the SN expands the competition 
between fading synchrotron emission from the emitting 
shell and decreasing free-free optical depth gives a char- 
acteristic light-curve with a relatively fast rise followed 
by a gradual decline, with the peak flux density occurring 
later at longer wavelengths. Those of our sources which 
have had rapid rises in flux density at 18 cm and stable 
or decreasing flux densities at shorter wavelengths are 
almost certainly radio SNe. The expected radio spectra 
of such sources is a power law with a sharp cutoff to- 
ward long wavelengths. While foreground ionized ISM 
can also cause low-frequency turnovers in SNR sources 
these turnovers are expected only at frequencies < 2 GHz 
(|Parra et al.l I2007D , hence sources with turnovers at a 
higher frequency are most likely SNe. Based on the 
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TABLE 3 

Source Resolution Probabilities and Size Estimates 



Source 


2 cm 


3.6 cm 




Radius 2 cm (mas) 






Radius 3.6 cm (mas) 




Radius 


Diameter 


Source ID 




(%) 


(%) 


Mean 


Mm 


Max 


Mean 


Min 


Max 


(mas) 


(PC) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Wll 


95.2 


68.8 


0.258 


0.142 


0.362 






0.331 


0.26 


0.19 


SN 


W12 


88.4 


24.8 






0.274 






0.201 


<0.28 


<0.21 


T 


W15 


99.8 


98.5 


0.358 


0.310 


0.418 


0.409 


0.308 


0.502 


0.38 


0.28 


SNR 


W17 


96.3 


90.5 


0.324 


0.202 


0.444 


0.410 


0.207 


0.584 


0.37 


0.27 


SNR 


W18 


99.9 


98.1 


0.527 


0.421 


0.629 


0.486 


0.328 


0.616 


0.51 


0.38 


SNR 


W25 


33.0 


72.4 






0.138 






0.285 


<0.29 


<0.22 


T 


W33 


96.6 


97.9 


0.317 


0.205 


0.429 


0.472 


0.320 


0.596 


0.39 


0.29 


SNR 


W34 


99.9 


4.8 


0.444* 


0.372* 


0.506* 






0.100 






SN 


W39 


92.7 


69.3 


0.309 


0.138 


0.472 






0.501 


0.31 


0.23 


U 


W42 


98.2 


99.1 


0.305 


0.225 


0.372 


0.544 


0.415 


0.672 


0.42 


0.31 


SNR 


W55 


97.7 


43.4 


0.202 


0.145 


0.229 






0.209 


0.20 


0.15 


SN 


W56 


92.5 


87.0 


0.177 


0.100 


0.226 






0.327 


0.18 


0.13 


SN 


W58 


99.8 


10.9 


0.389* 


0.317* 


0.450* 






0.101 






SN 


W60 


64.3 


9.8 






0.234 






0.100 


<0.24 


<0.18 


SN 


E10 


83.1 


97.6 






0.204 


0.280 


0.224 


0.328 


0.28 


0.21 


U 


Ell 


99.5 


97.4 


0.482 


0.367 


0.585 


0.500 


0.327 


0.656 


0.49 


0.36 


SNR 


E14 


99.8 


95.8 


0.460 


0.390 


0.519 


0.363 


0.227 


0.475 


0.41 


0.30 


SN 



Note. — Col. (1): Source name from Lonsdale ct al. (2006) for all sources except W55 and W56 which are from [Parra et al.l {2007). Sources W58 
and W60 arc newly detected sources. Cols. (2) and (3): probability (P) of source being resolved at 2 cm and 3.6 cm, respectively, using the method 
described in Section 13.2.11 We classify a source as provisionally resolved at that frequency in cases P > 90%. Cols. (4)-(6): the mean, minimum, 
and maximum outer shell radii at 2 cm (experiment GC028B) as calculated by the algorithm described in Scction |3.2.2| In cases where a source is 
unresolved only an upper limit is given, i.e., a max value. Cols. (7)— (9): the mean, minimum, and maximum outer shell radii at 3.6 cm (experiment 
GC02SA) as calculated by the algorithm described in Section 13.2.21 In cases where a source is unresolved only an upper limit is given, i.e., a max 
value. Col. (10): best estimate of outer shell radius. In cases where the source is resolved at both bands this is calculated by averaging the values of 
Columns 4 + 7 in mas and the result is printed in bold font. Where a source is resolved only at one band the radius is given as the mean value of 
the resolved band. Where a source is unresolved at both bands, the radius is given as an upper limit taken as the bigger of the two max values. Col. 
(11): outer shell diameter in parsccs. Col. (12): source identification; SN — supernova, SNR — supernova remnant, T — transition candidate, U — 
unclassified (see the Appendix). The * beside the 2 cm size estimates for W34 and W5S indicate that we consider these unreliable (sec Scction |3.2.3H 
and we therefore do not list a size for these sources in Columns 10 and 11. 



above light-curve and spectrum criteria in the Appendix 
we provisionally classify Wll, W34, W55, W56, W58, 
W60, and E14 as radio SNe. 

Transition objects. When the SN blast wave reaches the 
boundary of the wind-blown bubble it starts to interact 
with the constant density ISM. This boundary is deter- 
mined by the balance between wind ram pressure and 
ISM static pressure. As the shock wave propagates out- 
ward dense shocked ISM gas is accreted, relativistic par- 
ticles are e fficiently accelerated and the magnetic field 
is amplified (jBerezhko fc Volkll200l causing the source 
to brighten simultaneously at all radio frequencies. This 
transition phase lasts approximately until the swept up 
mass equals the ejecta mass and the source enters its 
SNR Sedov phase. In our data W12 and W25 have ris- 
ing light-curves at multiple wavelengths and are therefore 
candidate objects of this type. 

A similar increase in flux density at all radio wave- 
lengths has been obs erved in SN1987A in the Large Mag- 
ellanic Cloud (LMC; iZanardo et al.|[20Toh . Even though 
similar to the radio brightening expected in sources in 
which the blast wave starts to interact with the dense 
ISM, in the case of SN1987A it is thought that the SN 
blast wave has started to interact with a dense ring of gas 
emitted by the progenitor star 20, 000 years before the 
explosion (Bu rrows et al.l fl995). Hydrodynamic mod- 
eling suggests that the dense ring was emitted durin g 
the merger of two stars (jMorris fc Podsiadlowskill2007l ). 
It is possible that the potential "transition sources" in 
Arp220 also arise from a similar mechanism. Although 
such merger events must be rare in most galaxies they 
may be more common in the dense stellar environment 



of Arp220. 

Supernova remnants. When the shock created in an 
SN explosion has swept up an ISM mass equal to the 
ejecta mass the source enters the Sedov phase of SNR 
evolution. In the early part of this phase the total en- 
ergy in relativistic partic les stays approximately constant 
(jBerezhko fe ¥5113 72004). In models where the internal 
shell magneti c field comes from shock compressed exter- 
nal ISM field (van der LaaiJ[l96llThompson et al.ll2009ft 
the internal magnetic field is close to constant as the 
source expands and hence the radio luminosity is also 
constant. In contrast for m odels in which the ma gnetic 
field is internally amplified (|Berezhko fc Volkll2004[) , the 
magnetic energy density decreases and so does the radio 
luminosity. A constant or falling flux density at all wave- 
lengths is an expected signature of a radio SNR. In addi- 
tion, since the SNR shell is propagating through mainly 
neutral ISM local free-free absorption effects should be 
small and the radio spectrum close to a power-law. Al- 
though free-free absorption from foreground ionized ISM 
can occu r, this is likely to happen only at frequencies 
< 2 GHz (jParra et al.ll2007h . In the Appendix we classify 
sources W15, W17, W18, W33, W42, and Ell as likely 
SNRs based on the above expected properties. Most of 
these sources show relatively flat (a < 0.5) spectra from 
18 cm to 2 cm. Such spectral indices are unusual but not 
unprecedented for SNRs. Since all sources discussed in 
this paper must be detected at 3.6 cm and 2 cm, selection 
effects will bias spectral indices to small values. 

Active galactic nuclei (AGNs). It has been argued 
(jDownes fc Eck art 2007) that a supermassive black hole 
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Fig. 3. — Plot of estimated source outer s hell diameters at 2 cm and 3.6 cm. The estimated source sizes and error bars at each wavele ngth 
are assigned as described in Section I3.2.2I Symbol color indicates the results of source resolution testing described in Section I3.2.1I Red 
indicates that the source is resolved at both 2 and 3.6 cm defined as 99% confidence rejection of the null hypothesis that the source is a 
point. Green means that the source is provisionally resolved at 2 cm only and cyan that the source is provisionally resolved at 3.6 cm only. 
White represents sources that are unresolved at both bands. 



and an AGN are present within the western nucleus of 
Arp220. Radio emission from an AGN is likely present at 
some levels and so is potentially detectable. We would 
expect such an AGN associated compact source to be 
long-lived, but perhaps randomly variable in intensity. 
If as expected the radio emission arises within a jet then 
the spectrum should be flat (being generated from the 
superposition of synchrotron self-absorbed components 
peaking at different frequencies along the jet). Inspec- 
tion of Figure Q] shows that the 3.6 cm images of the 
flat spectrum sources W33 and W42 could be consistent 
with a jet morphology. Both sources are long-lived and 
the first shows strong 18 cm variability. More observa- 
tions, searching for structure and short-term variability, 
are needed to confirm the presence or absence of an AGN 
in Arp220. 

Unclassified sources. Two of the sources, described in de- 
tail in the Appendix, namely, W39 and E10 have proved 
impossible to classify because they show declining lumi- 
nosity at long wavelengths and increasing luminosity at 
short wavelengths, a behavior not predicted for any of 
the classes described above. More data are required on 
these sources. 

4.2. Flat Spectrum SNRs 

It was noted in Section 14.11 that the sources identified 
as SNRs mostly have flatter spectra than expected in 
standard models (with a — 0.5 — 0.7). Three alterna- 



tive models could explain the origin of such flat spec- 
tra with spectral ind i ces a < 0.5. In one such model 
Schlickcis er fc Fuerstl (|1989f ) argue that Fermi accelera- 
tion in strongly magnetized plasma flattens the electron 
injection spectrum (7 < 2) which consequently flattens 
the synchrotron emission spectrum giving a < 0.5. In an- 
other model the intrinsic spectral index is a ~ 0.5 in the 
SNR shell but the overall spectrum appears flatter due to 
spatially variable free-free absorption. The integration 
over the whole SNR of local spectra with turnovers due 
to free-free absorption happening at different frequencies 
produces a flatter global spectrum. In the last model the 
presence of a plerion-like component with a flat or in- 
verted spectrum in the center of the SNR results in an 
overall radio spectrum with spectra l index a < 0.5. Re- 
cent V LBI observations of SN1986J ()Bietenholz fc Bartell 
2008) show the emergence of a new radio component in 
the center of the expanding radio shell. The new compo- 
nent shows an inverted radio spectrum contrasting with 
the power law plus free-free absorption turnover spec- 
trum of the shell; the net result of this is to flatten the 
integrated spectrum. It is interesting that within Arp220 
the spectrum of source W33 might be "double humped" 
(see Figure [4]) showing a peak at ^2 GHz and another 
one at ^8 GHz. This is consistent with SN1986J like ob- 
jects. Adding further to this interpretation, it is observed 
that both W33 and W42 have significantly smaller size 
at 2 cm than at 3.6 cm (see Section [3.2.3)) which would 
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Fig. 4. — Spectra for the sources detected at 2 cm and 3.6 cm. Blue points are from measurements taken in time period 2006.02-2006.43 
and red points in time period 2006.91-2008.44. Error bars are plotted at ±lcr. In case of non detections upper limits are shown as triangles. 
At 18 cm (1.7 GH z), 13 cm (2.3 GH z), and 3.6 cm (8.4 GHz) upper limits are shown at upper limits of 6cr, 3.8<x, and 4.4tr respectively 
(see Section 3.1 of Parra ct al. 2007). The blue data points at 13 cm and at 3.6 cm were taken from simultaneous observations in BP129. 
The blue points at 6 cm (5 GHz) are taken from new reduction of BP129 observations. The new 6 cm, 3.6 cm and 2 cm (15 GHz) points 
are plotted in red. The 18 cm data are from unp ublished experiment GD21A (epoch 2006.43). The spectra are grouped according to 
the classification scheme discussed in Section 14.11 where SN = supernova, T = transition object, SNR = supernova remnant, and U = 
unclassified. For unresolved sources we adopt the measured peak brightness per beam as the best estimate of the total flux density to plot 
while for resolved sources we plot the integrated flux density within a tight box surrounding the source. 



be consistent with a compact central source being more 
dominant at higher frequency. 

4.3. Comparison of SNe and SNRs Sizes 

The top panel of Figure [5] plots the detected sources 
in the luminosity-diameter plane. This figure shows a 
clear difference in size between sources classified as SNRs 
(blue symbols) and SNe (red symbols). All six of the de- 
tected SNRs are resolved with diameter > 0.27 pc while 
all the SNe (except E14 which has size 0.30 pc) have 
sizes < 0.2 pc. In normal galac tic disks it has been 
argued (jChomiuk &i Wilcoti 120091 ) that the most com- 
pact/luminous SNR observed is determined by the age 
of the youngest SNR expected given the star formation 
rate. This mechanism cannot however explain the min- 
imum size of SNRs in Arp220 because given the high 
predicted SN rate (4 ± 2 yr _1 ) the above model predicts 
a much smaller limit than is observed. A physical min- 
imum size is in fact expected because SNRs should first 
"switch on" and become luminous just as they enter their 
Sedov phase, which occurs when the swept up ISM mass 
equals the ejecta mass. The radius at which a source 



enters the Sedov phase there fore depends on the ISM 
density and the ejecta mass. IScoville et al.l (|1997f ) esti- 
mate in the western nucleus a mean molecular density 
of 1.5 x 10 4 cm~ 3 ; more recently TSakamoto et al.l (2008) 
estimate from gas dynamics a mean total mass density of 
2 x 10 3 Mi 



1q pc 
sity < 4xl0 4 cm 



3 implying a gas molecular number den- 
3 . Adopting a molecular ISM density of 
10 4 cm -3 a minimum SNR diameter of 0.3 pc is expected, 
close to that observed. This cutoff size for SNRs goes 
only as the one third power of the assumed ejecta mass 
and external density and so is only weakly dependent on 
these quantities. Despite this the sharpness of the ob- 
served minimum cutoff in size suggests that the observed 
Arp220 SNRs are not embedded in densities > 10 5 cm~ 3 . 
This is either because such high density gas has low vol- 
ume filling factor in Ar p220 or because SN Rs become 
inefficient radio sources (| Wheeler et al.|[l980l ) when em- 
bedded in such high-density regions. 

In addition to a minimum size for SNRs we expect a 
maximum size for SNe set by the size of their wind-blown 
bubble. The predicted sizes of such bubbles are a func- 
tion of the external ISM pressure, progenitor mass-loss 
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rate and wind velocity (see Equation 2 of iParra et al.l 
120071) . The I SM pressure in Arp 220 is estimated to be 
10 7 K cm -3 (jDopita et al.l 120051 ). The wind mass- loss 
rate and speed are ex pected to vary significantly between 
different progenitors. I Weiler et al J (|2002f ) estimates how- 
ever that for the most powerful radio SNe (mostly Type 
Iln) mass-loss rates are 10~ 4 M yr _1 and velocities 
10 km s _1 . Adopting these parameters gives an esti- 
mated wind-blown bubble diameter of 0.4 pc consistent 
with the upper limits on SNe sizes and comparable in 
size to our SNR diameters. It should be noted however 
that there is considerable uncertainty in the expected 
wind-blown bubble size since only the ratio of mass-loss 
rate and wind velocity is constrained by VLBI SN obser- 
vations; and for a fixed ratio of these two quantities the 
bubble size is linearly dependent on the assumed wind ve- 
locity, a quantity which is poorly known. The fact that 
most of our observed SNe have size limits significantly 
smaller than 0.4 pc suggests that wind speeds could be 
5 km s^ 1 or less. 

In Section 14.11 we discuss the possibility of SN/SNR 
"transition objects" where interaction between the blast 
wave and the ISM has started but the amount of mass 
accreted is still less than the ejecta mass. Such objects 
would be expected to be intermediate in size between 
SNe and SNRs. The two transition object candidates we 
have, W12 and W25, have estimated diameters < 0.21 pc 
and < 0.22 pc, respectively, which is consistent with a 
progenitor with a somewhat lower mass- loss rate/wind 
velocity or a slightly higher ISM pressure than assumed 
above. 

It is notable that all of the SN candidates except one lie 
close together in the luminosity-size plane (see Figure [5]). 
All these positions would be reached by the evolutionary 
tracks of SNe more radio luminous by approximately a 
factor of two than the well-studied Type Iln SN1986J 
(shown by the thick black line marked 86 J in Figured]). 
One SN classified source (E14) is however exceptional 
since it lies at a position in the luminosity-size plane 
more characteristic of SNRs. Such an intense SN source 
which is still highly luminous at late times when it has a 
large diameter (and hence likely peaked at an even higher 
luminosity), would not however be unprecedented. The 
radio supernova SN1982aa and gamma ray bursts/SNe 
SN1998bw and SN2003w had 6 cm p eak luminosities 
nearly 10 times larger than SN1986J (|Chevalierl [20061) 
and their evolutionary tracks would be consistent with 
the position of E14 in the luminosity-size plane. 

4.4. Source Expansion Velocities 

Supernovae. The majority of the SN identified sources 
have insufficient data to estimate their explosion dates 
and hence set limits on their expansion velocities. Three 
sources E14, Wll, and W34 have however been detected 
at several 18 cm epochs showing close to linearly increas- 
ing flux densities with time (Batejat et al. 2012, in 
preparation). Comparing to equivalent portions of the 
fitted 18 cm light-curves for well-sampled radio SNe of 
similar luminosities, it seems that reasonable estimates 
of explosion dates (within accuracies of one year or so) 
can be made by linearly extrapolating the 18 cm light- 
curves down to zero flux density and then subtracting 
a year. Applied to E14, Wll, and W34, this algorithm 



gives ages, at the time of the 2 cm and 3.6 cm observa- 
tions, of 7, 6, and 6 years respectively. E14 is resolved at 
both 2 cm and 3.6 cm and Wll is provisionally resolved 
at 2 cm, allowing to measure their sizes with error bars. 
The size measurement of W34 we consider unreliable and 
cannot be used to estimate an expansion velocity (see 
Section 13.2.31) . Combining sizes and age estimates for 
the other two SNe we find that the E14 expansion veloc- 
ity is in the range 13, 700 - 28, 000 km s" 1 (90% confi- 
dence) and the Wll expansion velocity is in the range 
7300 - 15, 200 km s -1 . For the latter object the range is 
consistent with the canonical value of 10, 000 km s _1 for 
a normal Type II SN. The velocity for E14 is somewhat 
larger than this but would not be exce ptional. For in- 
stanc e VLBI observations of SN1993J (iMarcaide et alJ 
Il997f ) show it having a radio shell expansion velocity 
of 15, 000 km s _1 at early times. Also iBrunthaler et al.1 
(2010) have recently measured an expansion velocity of 
- 23,000 kms" 1 for SN2008iz in M82. 

Supernovae remnants. Of the six detected SNRs with 
measured sizes, four (W17, W18, W33, and W42) have 
been known since the original discovery observations of 
ISmith et all (|1998f ) made in 1994. The remaining two 
sources, W15 and Ell, have stable 18 cm light-curves and 
spectra consistent with an SNR origin (see the Appendix) 
but with luminosi ties below the detection limit of the 
Smit h~et~aT] (|1998l ) observations; these sources were very 
likely also present in 1994 but were not detectable. It 
seems highly probable that all the SNR identified sources 
are at least 12 years old . In fact 18 cm ligh t-curve mon- 
itoring (see Figure 3 in lRovilos et al.ll2005| ) shows their 
light-curves declining relativel y slowly so they are li kely 
considerably older than this. iLonsdale eFail (I2006T) es- 
timate ages of several decades. Our recent reanalysis of 
all the 18 cm monitoring data (Batejat et al. 2012, in 
preparation) agrees with this conclusion. Taking a rough 
estimate of their ages as 20 years and taking the mean 
diameter for this group of sources (0.31 pc) gives average 
expansion speeds since explosion of 5000 km s _1 . Such 
speeds are consistent with those expected for SNRs just 
entering the Sedov phase when the swept up ISM mass 
equals that of the ejecta so that strong deceleration is 
occurring. 

4.5. Source Magnetic Fields and Energetics 

Given estimates of source sizes and synchrotron emis- 
sion flux densities it is possible to work out minimum 
energies and equipartition magnetic fields. Most mod- 
els of SNR evolution predict however that SNRs are far 
from equipar tition and that energies are particle dom- 
inated. The iBerezhko fc Vo lk (2004) model assumes a 
1% ratio of magnetic field to relativistic particle energy 
density. For the seven sources which are well resolved 
at two bands (comprising six SNR candidates and one 
SN candidate) we give in Table @] estimates of magnetic 
fields and energies calculated assuming the revised mag- 
n etic field - relativistic p article equipartition expression 
of lBeck fc Krausel (|2005[ ) adjusted to a 1% magnetic field 
to particle energy density ratio assumption. We assume 
a spherical shell geometry with a ratio of inner to outer 
radius of 0.8, an outer shell diameter given in Table [3] 
and an internal volume filling factor for radio emission of 
10%. We also assume a ratio of proton to electron num- 
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TABLE 4 

Magnetic Field and Energy Estimates 



Source 


a 


Diameter 


B Rev 


E Part 


E Part, DW 






(pc) 


(mG) 


(xlO 49 erg) 


(xlO 49 erg) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


W15 


0.72 


0.28 


23.4 


3.3 


3.5 


W17 


0.23 


0.27 


46.7 


12.0 


3.1 


W18 


0.62 


0.38 


19.9 


6.0 


1.8 


W33 




0.29 


18.0 


2.2 


3.9 


W42 


0.24 


0.31 


49.7 


20.4 


1.4 


Ell 




0.36 


15.4 


3.1 


5.2 


E14 


0.72 


0.30 


20.9 


3.3 


1.0 



N ote. — Col. (1) : source name. Col. (2): spectral index taken 
from Parra et al.l (2007) . The revised cquipartition formula is not valid 
for spectral indices < 0.5. In cases where a source has a < 0.5, i.e., 
sources W17 and W42, a is taken as 0.51. In cases where no spectral 
index information is available, i.e., sources W33 and Ell, a is taken as 
0.7. Col. (3): source diameter taken from Table [3] Col . (4): magnetic 
field d erived using the revised equipartition formula of Beck & Krausc 
(2005) adjusted to a magnetic field to particle energy ratio of 0.01. K, 
the ratio of number densities of cosmic-ray protons and electrons per 
particle energy interval within the energy range traced by the observed 
synchrotron emission is taken as 100. A shell inner to outer radius of 0.8 
and a synchrotron filling factor of 10% arc assumed. Col. (5): associated 
total energy in rclativistic particles. Col. (6): total energy in rclativistic 
particles estimated assuming energ y density bala nce with ram pressure, 
expansion velocity from Drainc &c Woods (1991), a molecular number 
density of 10 4 cm -3 , and a synchrotron filling factor of 10%. 



be r density at fix e d ene rgy K = 100 which according 
to IBeck fc Krausel (pOOl is valid for young SNRs. The 
results give magnetic fields in the range ~ 15-50 mG 
and total energies in relativistic particles between 2% 
and 20% of the expected 10 51 erg kinetic energy of a 
typical SN. This is consistent with estimates given by 
lLacki et al.l (|2010| ) for total particle energy fractions of 
relativistic particles needed to explain the FIR-radio cor- 
relation in compact star bursts via calorim e ter m odels. 

Consistency with the iBerezhko fc Volkl (|2004D model 
can be checked by comparing the particle energy 
in an SNR with the prediction that it should be 
# P art = /0iSMV 2 Vsheii, where pism is the average density 
of the ISM, v is the SNR expansion velocity and V^hcii is 
the synchrotron emitting volume of the SNR. This calcu- 
lation assumes that ram pressure balances shell internal 
pressure which in turn is dominated by relativistic par- 
ticles. Consistent with the latter assumption there is 
increasing observational evidence for a large/dominant 
fractio nal relativistic particle e n ergy density in youn g 
SNRs (jPatnaude fc Fesenl [200l IBerezhko e'tldl [2009h . 
In Column 6 of Table 21 we give the estimated rela- 
tivistic particle density given by the above ram pres- 
sure balance formula. A source-emitting volume V^heii 
was calculated from the measured source outer diame- 
ter assuming a shell with ratio of inner to outer radius 
of 0.8 with synchrotron filling factor of 10%. The ram 
pressure pismv 2 was ca lculated from velocities given by 
IDraine fc Wo"odsl ([l99l for SNRs of the size observed in 
a medium of density 10 4 cm -3 ; note however that this 
estimate is only weakly dependant on the assumed ISM 
density since in the Sedov phase the predicted expansion 
velocity scales as p~ 5 for fixed diameter. As expected 
during the Sedov phase the estimated particle energies 
in Column 6 of Table 2] are close to a constant. We find 
that these estimated relativistic particle energies are in 
good agreement with those in Column 5 derived from the 



measured synchrotron flux densities assuming a 1% ratio 
of magnetic to particle energy density. 

Estimates of source magnetic fields are given in Col- 
umn 4 of Table [4] as derived from the measured source 
radio luminosities and sizes assuming a 1% magnetic 
to relativistic particle energy density ratio and the 
source geometry and synchrotron volum e filling factor 
described above. I Thompson et al.1 (|2009D estimate some- 
what smaller field values, averaging 9 mG for the sources 
listed in Table [4] These published estimates were made 
without the benefit of knowing the source sizes by assum- 
ing that 1% of the total kinetic energy of 10 51 erg goes 
into relativistic electrons and finding the magnetic field 
that gives the observed source radio synchrotron lumi- 
nosity (vL v ). It was argued that these relatively modest 
fields could be produced by shock compression of ISM 
fields of a few mG; fields sufficiently strong to explain the 
IR-radio correlation via calorimeter models, even in the 
presence of inver se Compton energy lo sses off ambient 
starlight photons (j Thompson et a l. 2009). This compres- 
sion model predicts however that SNRs have constant ra- 
dio intensities until the electron synchrotron energy loss 
time (> 200 years at 5 GHz for the sources in Arp220). 
Most observable sources would then be older than 100 
years old and therefore be larger than 0.7 pc in diamete r 
in a 10 4 cm~ 3 density medium (jDraine fc Woods! fl99l . 
a factor two greater than their measured sizes. An ad- 
ditional consequence of such a model, given the high SN 
rate in Arp220, is that many more SNR sources than the 
six candidates presented in this paper would be expected 
to be detected (order of a few hundred) . The low number 
of observed SNRs in Arp220 combined with their small 
sizes argue against models in which SNR magnetic fields 
are generated from compressed ISM magnetic fields. 



4.6. The SNR Luminosity- Size Relation 
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A correlation between the radio surface brightness (£) 
of SNRs and the ir diameters (D) has long been claimed 
(Shklovskii 1960). Such a correlation can alternatively 
be cast as one between radio luminosity (L) and di- 
ameter (£>), a formulation which removes the implicit 
D ~ 2 correlation automatically induced by plotting S. 
Both the reality and slope of the correlation has been 
a controversial topic given the likeli hood of stron g se- 
lection effects, both in our Galaxy (IGreenl [2001 and 
(to a lesser extent) in external galaxies (|Urosevic et al.l 
2005). It has also been claimed that if the £-D and L- 
D correlations exist, they are seconda ry correlations due 
to correlations with the I SM density (jBerkhuiisenl 1 1 9861 : 
iBandiera fe Petrukl 120101). Finally, it has been claimed 
(jArbutina fe Urosevid 20051 ) that different correlations 
exist for the luminous SNRs in galactic dense clouds and 
most observed extragalactic SNRs compared to galactic 
SNRs in low density regions. 

We can test the claimed L-D correlation for powerful 
SNRs by comparing the diameters and luminosities of the 
SNRs in Arp220 with those found in M82 and the LMC 
(see bottom panel of Figure [5]). The data for 45 SNRs in 
M82 used in Figure[5]are taken from lHuang et all (|1994f ) 
who use the FWHM of a Gaussian fit as an estimate 
of source diameter. To make these estimates compati- 
ble with our Arp220 estimates we have converted this to 
a dispersion and then to the diameter of a shell model 
with ratio of inner to outer radiu s 0.8. The data for the 
LMC are taken from iMills et al.l (| 19841 ). These authors 
obtained sizes by taking two models, one comprising a 
thin circular ring and another a uniform circular disk 
and fitting to the half power response along the major 
and minor axes. An average of the two values is chosen 
and is expected to give a reasonable approximation to 
the shell size as defined for Arp220. 

Figure [5] clearly shows that the data for Arp220 are 
to first order consistent with an extension of the L-D 
correlation to very small size SNRs. More quantitatively 
an unweighted least-squares fit to all SNR points gives 
the best-fitting relation L oc D~ 19 (drawn in blue in 
Figure [5]) . This gradient is however mostly determined 
by the M82 and LMC data because of the larger number 
of SNRs plotted for these galaxies. 

Models f or radio SNR evolu tion s uch as those pre- 
sented by iHuang et all ([ 19941 ) and iBerezhko fc Volkl 
(2004) posit constant ratios between relativistic par- 
ticles and field energy densities and roughly con- 
stant energy in relativistic particles during the Scdov 
phase. Individual SNRs evolve along the luminosity- 
diameter correlation during their Sedov phase on tracks 
which are independent of the external density nisM 
(|Berezhko fc Volkl 120041 ) PI Assuming an optically thin 
radio synchrotron spectral index of 0.5, these models pre- 
dict L oc D~ 9 / 4 E 7 / 2 , where E is the SN kinetic energy. 
Assuming that all SNe have close to the same energy 
then we expect L oc £)~ 9 / 4 = D~ 2 25 . Assuming this de- 

Despite this direct independence of L on riigM these models 
can explain that the most luminous SNRs are found in host nuclei 
with high densities such as Arp220 because (1) such nuclei have 
the highest star formation rates and h ence the youngest, smallest 
diameter, and so most luminous SNRs (Chomiuk & Wilcots 2009) 
and because (2) there is a minimum "switch on" size for Sedov 
phase SNRs, which decreases as density increases, increasing the 
maximum SNR luminosity. 



pendence and fitting the constant of proportionality we 
obtain the green line in Figure [5j This theoretical corre- 
lation fits the observed Arp220 data very well and also, 
given their large internal scatter, is consistent with the 
M82 and LMC data. 

Some competing models argue that radio SNR lumi- 
nosity is instead mostly determined by ISM density with 
source size being of secondary importance. In o ne such 
model ()van der Laa n 1962; Tho mpson et al.l2009f ) the in- 
ternal magnetic fields that control SNR luminosity are 
compressed ISM fields which themselves increase with 
ambient density (£?ism oc nf gM ) giving rise to a L-nisM 
relation. A secondary L-D relation can however still 
occur because the diameter of typically observed SNRs 
also depends on ISM density. Specifically in the model of 
Ivan der Laanl (1962) the radio luminosity of an SNR at a 
given frequency stays roughly constant till a time equal 
to the electron synchrotron loss time, t syn after which 
it rapidly declines. Given that SNRs are decelerating 
during their Sedov and subsequent pressure-driven snow- 
plow phases it follows that most observed radio SNRs in 
a flux-limited sample will have sizes close to their sizes 
at the synchrotron loss time t syn . For most sources plot- 
ted in Figure O the predicted i syn is comparable or larger 
than the age and which they enter the snow-plow phase 

for which D oc t%n (jMcKee fc Ostrikerlll977D ; sub- 

stituting t syn oc -Bjsm anc ^ ftsM <x n \sM a -C~ n iSM re- 
lationship results. Combining this with the predicted 
L—nisM relationship noted earlier we find that we get an 
L-D relationship with power-law exponent of —2.25 if 
the power-law exponent linking ISM magnetic field and 
d ensity has a value a = .6. 

iThompson et all (120090 assuming the Ivan der Laanl 
(1962) SNR model have estimated -Bism in a num- 
ber of galaxies and have compared this to estimates 
of ISM surface density £ finding the empirical relation 
([Thompson et al.ll2009l Equation (6)) B oc S 55 . If one 
assumes a constant gas scale height then n oc S and 
-Bism oc n 055 ; similar to the value of the exponent re- 
quired above to explain the slope of the L-D relation. 
Despite this consistency there is no physical explana- 
tion of why the exponent value is .55rj In contrast 
the model of [Berezhko & V61k ( 2004]) directly predicts 
the slope of the observed L-D from the model's assumed 
physics. Additionally in Section 14.51 we gave arguments 
based on source energet ics, sizes and numbers that the 
IBerezhko fc Volkl (120041) model applie s in Arp220. The 
fact that the IBerezhko fc Volkl (|2004D model naturally 
explains the L — D relation between galaxies and that all 
galaxies fall on the same L-D relation argues that this 
is the only mechanism n eeded and that it also ap plies in 
M82 and the LMC. The IBerezhko fc Volkl (pOM ) model 
in normal galaxies is al so favored by the work on SNR 
luminosity functions bv lChomiuk fc Wilcotsl (|2009f ). 

A strong test which can diff erentiate between d ensity - 
dependent models and the IBerezhko fc Volkl (|2004f ) 
model would be to look for the time variations in lu- 
minosity predicted to occur only in the latter model as 
individual sources expand and move along the L-D cor- 

2 |Thompson et al. (2009) instead argue that physically 
Bism oc n might be expected based on equipartition between the 
magnetic field energy density and the energy density required to 
balance disk gravity 
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Fig. 5. — Illustration of the 3.6 cm flux densities and sizes of observed SNe and SNRs. The bottom panel plots sources in Arp220 and 
SNRs in M82 and the LMC. The top panel is a close up view of the dashed rectangular region shown in the bottom panel showing only 
Arp220 sources. Red is used for sources identified as SNe, blue for sources identified as SNRs, green for possible transition objects, and 
cyan for unclassified source s The SN tracks fo r SN1986J (Type Iln), SN1979C (Type IIL), and SN1980K (Type I I L) we re produced using 
the light-curve fits given in Wcilcr ct al. (2002) combined with the deceleration parameter from Bictcnholz ct al. (2002) for SN1986J and 
assuming free expansion at 10 4 km s _1 for both SN1980K and SN1979C. The square markers along each track indicate time evolution 
and are 1 year apart. Black crosses and dots show SNRs in M82 and the LMC. W34 and W58 are omitted as their size estimates are not 
considered dependable. In cases where a point is den oted by a triangle pointing left this signifies that the size estimate is an upper limit. 
In both panels the green line is the theoretical line of Bcrczhko & Volk (2004) of slope —9/4 whi le the blue line is th e best fit through the 
SNR luminosities and s izes for all three g alaxies. The dat a for the SNRs in M82 are taken from [Huang et aLl i 19941 ) and the data for the 
SNRs in the LMC from lMills et al.l lfl98l . see Section |4~6)| . 

galaxies. The compressed ISM magnetic field contri- 
bution must therefore be significantly less than our to- 
tal estimated internally generated field, i.e., /-Bism <C 
- Brnr, with / ra nging from 3 to 6 for young SNRs 
( Volk et al.ll2002f L Taking the median 5 S nr horn Ta- 
ble|4]we get an upper limit on Bism ranging from 3.5 mG 
to 7 mG in Arp220. This is consistent with the esti- 
mates of magnet ic fields from OH ma ser Zeeman split- 
ting (0.7-4.7 mGjRobishaw et al]|2008ft . For comparison 
iThompson et a l. (2009) estimate a minimum ISM mag- 
netic field of ~ 1 mG in order that sync hrotron emissivity 
domin ates over inverse Compton losses. IThompson et all 
( 2009) also estimate an upper limit of 20 mG based on the 
argument that the magnetic field is dynamically lower 
than gravity. 

5. CONCLUSIONS 
The main conclusions of this paper are as follows. 



relation. Given the small sizes of the SNR sources in 
Arp220 quite large flux variations are expected over rel- 
atively short times. For instance if our SNRs are expand- 
ing at 3 000 km s _1 just as they enter the Sedov phase, 
then for a diameter ~ 0.3 pc they increa s e in d iameter by 
2% yr -1 which from iBerezhko fc Volk! ([2 004) im plies a 
flux density decrease of 4.8% vr- 1 . IRovilos et al] ((2005T) . 
who analyzed five epochs of 18 cm data spanning 5.6 
years, were able to rule out factor of two variations over 
that period, as expected from SNe models, but not varia- 
tions of the amplitude predicted above. We are presently 
analyzing nine 18 cm data sets over a longer time period 
to see if we can detect the predicted flux density decline. 

4.7. ISM Magnetic Fields in Arp220 

In Sections 14.51 and 14.61 we argue that in SNRs in- 
ternally generated magnetic fields dominate over com- 
pressed ISM magnetic fields both in Arp220 and other 
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1. We have detected two new sources at both 2 cm 
and 3.6 cm wavelength in the western nucleus of Arp220, 
namely W58 and W60. 

2. We have resolved for the first time, 11 of the 17 de- 
tected sources at 2 cm, 8 at 3.6 cm, and 7 at both 2 cm 
and 3.6 cm. 

3. We confirm the claim of iParra et all ()2007[ ) that 
the compact radio sources in Arp220 mostly comprise 
a mixed population of SNe and SNRs. Two sources are 
candidate SN/SNR transition objects. A few sources re- 
main difficult to classify and may be AGN components. 
Ongoing VLBI monitoring at 6 cm should shed further 
light on the nature of these sources. 

4. The sources resolved at both wavelengths (all but one 
of which are SNRs) have diameters in the range 0.27- 
0.38 pc with mean 0.31 pc. In comparison the SNe are, 
except in one case (E14), all unresolved. The observed 
size boundary between SNe and SNRs is consistent with 
an ISM density of - 10 4 cm" 3 . 

5. Combining source sizes with source ages enables us 
to calculate upper limits for source expansion velocities 
averaged over their lifetime. We find expansion veloc- 
ities < 30,000 km s" 1 for the SNe Wll and E14 and 
~ 5000 km s- 1 for the SNRs. 

6. We argue that magnetic fields in the SNRs o f 
Arp220 are internally generated (Bcrezhko & Volk 2004) 
and not dominated by compressed ISM magnetic fields 
(|van der Laanl I1962D . This interpretation is supported 
by the fact that the particle energy density estimated 
from synchrotron fluxes and sizes, assuming a 1% ratio 
of magnetic field to relativistic particle energy density 
ratio, is equa l to the ram pressure ene rgy density, just as 
predicted by iBerezhko fc Volkl (|2004f) . Furthermore the 
relativistic particle total energies are consist ent with val- 
ues r equired to fit the FIR-radio correlation (jLacki et al.l 
120101 ). In contrast the compressed field model is incon- 
sistent with the low number of observed SNRs and their 
small measured diameters. 

7. The observed SNR radio luminosity as a function of 
diameter for sources in Arp220, M82 and the LMC is in 
good agreement with the predicted relation L v oc £)~ 2 - 25 
as derived ass uming internally generat ed magnetic fields 
in SNR shells (IBerezhko fe Volkl[200l . 

8. Based on our conclusion that in SNRs internal mag- 
netic fields dominate over compressed ISM magnetic 
fields we estimate an upper limit of 7 mG for the ISM 
magnetic field in Arp220. This is consistent with other 
estimates and limits. 
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APPENDIX 

This appendix assigns source type classifications 
(mostly SN or SNR) to the high frequency detected 
sources in Arp220. In Section 14.11 we discuss the 
expected temporal and spectral properties of each of the 
possible source classes. The most critical information 
for classification is the variability properties of each 
source. Concentrating on recent high sensitivity data we 
find that source variability properties fall naturally into 
four different groups which we term "rapidly rising" , 
"possibly rising", "stable", and "possibly declining". 
The term "rapidly rising" is used in the individual 
source descriptions below at 18 cm and 6 cm for a source 
whose flux density has risen more than 3<r and 30%, 
between two high sensitivity epochs over ~ 2.5 years 
(comparing the two 6 cm epochs in Table [2] and the 
2003.85 and 2006.44 epochs at 18 cm). The same term is 
used at 3.6 cm for a source whose flux density has risen 
by more than 3<r and 15% over the 0.89 years between 
the two epochs. The term "possibly rising" is used for 
a source whose flux density has risen by more than la 
and the term "possibly declining" is used for a source 
whose flux density has declined by more than la. The 
term "stable" is used for any variation in flux density 
lower than la. 

Wll - First detected at 18 cm in epoch 2003.85 this 
source is rapidly rising at 18 cm. It is also rapidly rising 
at 6 cm and stable at 3.6 cm. Its spectrum peaks around 
3 GHz. We classify Wll as an SN. 

W12 - First detected at 18 cm in epoch 2002.88 this 
source is rapidly rising at 18 cm and 6 cm. A similar 
rise rate at 3.6 cm is consistent with the data. The radio 
spectrum shows a straight power-law spectrum from 
18 cm to 3.6 cm. We classify this source as a candidate 
SN/SNR transition object. 

W15 - First detected at 18 cm in epoch 2002.88. 
IParra et al.1 (|2007l ) classified this source as ambigu- 
ous because it was detected in the 18 c m observations 
GD17 A and GD17B but at levels below the lRovilos et al.l 
(2005) sensitivity limit, meaning that it could be either 
a new source or a long-lived stable source not detected 
earlier because of sensitivity limitations. Subsequent 
data, showing a stable 18 cm light-curve, are more 
consistent with the latter interpretation. However, a 
possible rise at both 6 cm (16%) and 3.6 cm (18%) may 
be in contradiction with an SNR origin. Despite this, 
considering potential ~ 10% calibration/reconstruction 
uncertainties, the large size of this object (0.28 pc re- 
solved at both 6 cm and 3.6 cm bands) and its spectrum 
peaking around 3 GHz, we provisionally classify W15 as 
an SNR. 

W17 - A long-lived source discovered by iSmith et al.1 
(|1998[ ). Its long-term 18 cm light-curve is stable and its 
6 cm and 3.6 cm light-curves are possibly declining. Its 
spectral shape is quite flat but could be interpreted as 
being double humped with one peak at ~ 1.5 GHz and 
another at ~ 5 GHz. We classify W17 as an SNR. 

W18 - A long-lived source first discovered by lSmith et al.l 
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(1998). Long term 18 cm flux density monitoring shows 
a steady decrease of 8% yr" 1 over 11.6 years. In more 
recent shorter wavelengths at 6 cm and 3.6 cm the 
source is classified, respectively, as possibly declining 
and possibly rising. W18's spectrum turns over around 
3 GHz. We provisionally classify W18 as an SNR. This 
classification is based primarily on its long-lived nature 
and declining light-curves at 18 cm and 6 cm. 

W25 - First observed at 18 cm in epoch 2002.88. It is 
rapidly rising at 18 cm, 6 cm and 3.6 cm. Its spectrum 
peaks around 3 GHz. Primarily based on its rapidly 
rising light-curves at all frequencies we consider this 
source to be a candidate SN/SNR transition object. 

W33 - A long-lived source first discovered in lSmith et al.l 
(1998). Its 18 cm and 6 cm light-curves are possibly 
declining over ~ 2.5 years while its 3.6 cm flux density 
has remained stable over 0.89 years. On average its 
18 cm flux density has decreased ~ 10% yr" 1 since 
monitoring began in 1994. The spectrum of W33 is 
quite flat with hints of two broa d peaks at ^2 GHz 
and ~ 8 GHz. iParra et al.l (|2007f ) suggested that W33 
is a possible AGN candidate. Inspection of the 3.6 cm 
image of W33 does not exclude the possibility of a jet 
morphology. Based on the behavior of its multifrequency 
light-curves we provisionally classify W33 as an SNR 
but an AGN origin is still possible. 

W34 - First discovered at 18 cm in epoch 2003.85. This 
was followed by a rapid rise at 18 cm epochs. W34's 
6 cm light-curve is also rapidly rising while the source 
is stable at 3.6 cm. W34's spectrum turns over around 

5 GHz. The light-curves and spectrum are consistent 
with it being an SN. 

W39 - A long lived source detected by ISmith et al.l 
(1998). Long term monitoring shows a steady decline at 
18 cm (~ 7% yr -1 ). From recent short wavelength data 
it is classified as stable at 6 cm and possibly rising at 
3.6 cm. The source's spectrum turns over at ~ 2 GHz. 
Given its unusual multi-frequency variability properties 
(declining at long wavelengths but stable or rising at 
short wavelengths) we refrain from speculation about 
the nature of W39 until more data are available. 

W42 - A long-lived source first detected by ISmith et al.l 
(1998). It has decreased by less than 8% in 18 cm flux 
density since 1994 while its recent 6 cm and 3.6 cm 
light-curves are stable. W42 is the most stable of the 
long lived sources. The spectrum is quite flat from 
18 cm to 3.6 cm. There is no sign of a low frequency 
turnover but there is a high-frequency turn down from 
3.6 cm to 2 cm. We conclude that W42 is likely to be a 
SNR. 

W55 - First detected at high frequency by IParra et al.l 
(2007) in epoch 2006.02. No 18 cm data are available. 
The source's 3.6 cm light-curve is possibly declining, the 

6 cm light-curve is rapidly declining and its spectrum 



peaks at around 8 GHz. We provisionally classify W55 
as an SN. 

W56 - This source was first discovered by IParra et al.l 
(|2007t) in epoch 2006.02. No 18 cm data are available. 
The source is rapidly rising at 6 cm and stable at 
3.6 cm. W56's spectrum turns over around 5 GHz. Con- 
sidering these properties we conclude that W56 is an SN. 

W58 - A newly detected source in the most recent 2 cm 
and 3.6 cm observations. At these two wavelengths it 
has flux densities of 226 and 322 /iJy, respectively. At 
this 3.6 cm flux density the source would have been less 
than a 4cr rms detection in BP129 (epoch 2006.02). New 
6 cm data from GC031A yield a flux density of 723 /iJy. 
Re-reduction of BP 129 6 cm data made it possible 
to clearly detect this source at epoch 200 6.02 with a 
flux density of 408 ^tJy (not detected by IParra et al.l 
12007ft allowing us to classify it as rapidly rising (77% 
increase in ~ 2.5 years). No 18 cm or 3.6 cm light- 
curve information is available. We classify W58 as an SN. 

W60 - A newly detected source in the most recent 
2 cm and 3.6 cm observations. The source was also 
detected in the latest 6 cm observation in experiment 
GC031A (epoch 2008.44) at a flux density of 1016 /LtJy. 
Re-reduction of BP129 6 cm data made it possible to 
clearly detect this source at epoch 2 006.02 with a flux 
density of 285 ^Jy (not detected by IParra et ail 12007ft 
allowing us to classify this source as rapidly rising at 
6 cm. We classify W60 as an SN. 

E10 - First detected at 18 cm in epoch 2002.88. It 
showed a decline in the next epoch and was undetectable 
in the last two 18 cm epochs. The source spectrum peaks 
at ~ 5 GHz. This observed fast low frequency decline 
is too rapid to be consistent with an SNR and suggests 
instead a rapidly evolving Type Ib/c SN observed after 
its peak at 18 cm. However in contradiction to this inter- 
pretation the source is rapidly rising at 6 cm and 3.6 cm. 
We leave E10 unclassified until further data are available. 

Ell - First detected in the high sensitivity 18 cm epoch 
2002.88 at a flux density below the detection limit of 
earlier epochs. Its 18 cm flux density has decreased by 
less than 4% in subsequent 18 cm monitoring spanning 
^3.5 years. It was not detected by IParra et al.l (|2007T ) 
at 6 cm and 3.6 cm. Re-reduction of BP129 6 cm data 
(epoch 2006.02) made it possible to clearly detect this 
source with a flux density of 429 ^Jy compared to 310 
/zJy in experiment GC031A (epoch 2008.44) allowing us 
to classify it as possibly declining at this wavelength. 
Based on its stable 18 cm light-curve and possibly 
declining 6 cm light-curve we classify Ell as an SNR. 

E14 - First detected at 18 cm in epoch 2002.88. In 
subsequent monitoring it shows a rapidly rising 18 cm 
light-curve over ~ 2.5 years. The source is classified as 
possibly declining at 6 cm and stable at 3.6 cm. The 
spectrum shows a turnover frequency at 2 GHz. We 
classify E14 as an SN. 
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